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ABSTRACT 


A  literature  survey  was  made  of  the  experimental  values  of  the  specific 
heat  of  graphite.  Most  of  the  measurements  from  20*K  to  3800#K  were  re¬ 
analyzed  and  a  new  average  curve  is  given  for  each  experimental  run.  Some 
of  the  new  curves  are  significantly  different  from  the  original  curves.  Esti¬ 
mates  are  given  of  the  Cp-Cv  term,  the  electrcnic  specific  heat,  and  the 
specific  heat  of  the  lattice  vibrations.  The  lattice  specific  heat  was  approxi¬ 
mated  by  various  combinations  of  Einstein  and  one-  and  two-dimensional 
Debye  functions.  The  characteristic  temperatures  were  selected  by  a  least- 
squares  curve-fitting  procedure. 

This  report  has  been  reviewed  and  is  approved. 
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1.  INTRODUCTION 


This  report  gives  the  results  of  a  literature  survey  of  the  experimental 
values  of  the  specific  heat  of  graphite  and  the  results  of  attempts  to  approxi¬ 
mate  the  lattice  specific  heat  by  relatively  simple  analytic  functions.  The 
survey  covers  most  of  the  experimental  values  of  the  specific  heat  at  constant 
pressure  Cp  from  20  to  3800° K.  The  lattice  specific  heat  Cj  is  estimated 
by  subtracting  from  the  average  experimental  values  of  Cp  an  estimate  of  the 
thermodynamic  formula  for  (Cp  -  Cv)  and  an  estimate  of  the  specific  heat  of 
the  electrons  Ce.  The  lattice  specific  heat  was  approximated  by  various 
combinations  cf  Einstein  and  one-  and  two-dimensional  Debye  functions.  The 
Debye  or  characteristic  temperatures  were  selected  by  a  least-squares 
curve-fitting  program  on  an  electronic  computer. 

Several  of  the  existing  literature  surveys  (,-s)  of  thermal  properties  in¬ 
clude  the  specific  heat  of  graphite.  The  present  study  differs  from  these  in 
two  important  respects.  First,  the  original  experimental  data  have  been 
reanalyzed,  whenever  possible,  and  an  average  specific  heat  curve  obtained 
for  each  sample  measured.  Second,  these  average  specific  heat  curves  for 
individual  samples  are  presented  in  tabular  and  graphical  forms  which  permit 
an  evaluation  of  the  experimental  uncertainty  in  the  average  specific  heat  for 
all  samples  and  which  show  to  what  extent  the  specific  heat  differs  for  different 
types  of  graphite.  The  original  experimental  data  were  reanalyzed  for  several 
reasons,  the  simplest  being  to  provide  a  check  against  numerical  errors  in 
data  reduction.  Much  of  the  early  literature  gave  the  mean  specific  heat  over 
large  temperature  intervals.  In  these  cases  the  true  specific  heat  at  each 
temperature  has  oeen  calculated  from  the  original  enthalpy  measurements. 
Some  of  the  recent  literature  report  specific  heat  curves  with  sharp  bends  or 
with  sections  which  decrease  with  increasing  temperature.  Such  behavior 
seems  unlikely  and  appears  to  the  present  author  to  be  due  to  the  manner  by 
which  the  specific  heat  was  obtained  from  the  enthalpy  measurements.  The 
total  specific  heat  function  is  the  sum  of  six  complicated  analytic  functions 
representing  the  contributions  of  the  six  branches  of  the  frequency  distri¬ 
bution  in  "Tc  space.  It  appears  to  be  impossible  to  accurately  represent  the 
specific  heat  function  over  a  several  hundred  degree  temperature  interval 
by  only  2  or  3  terms  of  a  power  series  expansion  in  the  temperature.  On  the 
other  hand,  when  more  terms  arc  used  and  the  coefficients  arc  determined  by 
a  least-squares  curve  fitting  to  the  experimental  enthalpy  data,  then  in  practice 
the  resulting  curve  tends  to  follow  small  drifts  in  the  data  which  result  in  the 
erroneous  features  of  the  specific  heat  curve  mentioned  above.  Consequently, 
it  appears  that  "best  by  eye"  curve  fitting  yields  all  the  accuracy  which  can  be 
obtained  in  most  cases,  and  this  method  has  been  used  exclusively  in  this  work. 

For  many  purposes  it  is  sufficient  to  have  the  empirical  specific  heat 
curve  in  tabular  form.  However,  for  certain  theoretical  studies  and  for  some 
computer  numerical  computations  it  is  moro  convenient  to  represent  the 
empirical  specific  heat  curve  by  an  analytic  function.  In  1923  Magnus'9*  used 
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a  sum  of  two  three-dimensional  Debye  functions  with  Debye  temperatures  of 
760  and  2280*K  to  approximate  the  specific  heat  curve.  The  analytical  and 
empirical  curves  crossed  at  several  temperatures  but  between  these  tempera¬ 
tures  the  difference  was  usually  much  larger  than  the  experimental  uncertainty 
in  the  empirical  curve.  Later  it  was  recognized  that,  because  of  the  weak 
binding  between  layers,  graphite  should  resemble  a  two-dimensional  crystal 
and  its  specific  heat  might  be  better  approximated  by  two-dimensional  Debye 
functions,  for  which  the  specific  heat  is  proportional  to  T2  at  low  temperatures. 
Still  later  Krumhansl  and  Brooks  (*•)  pointed  out  that,  as  the  temperature  ap¬ 
proaches  the  absolute  zero,  the  Tz  variation  must  change  to  a  T3  variation. 

The  trend  of  the  specific  heat  curve  in  this  transition  region  is  now  satisfactorily 
understood  from  the  results  of  a  number  of  theoretical  calculations'11) but  the 
actual  formulas  are  rather  complicated.  Since  the  objective  of  the  present 
work  is  to  get  tractable  formulas  to  approximate  the  specific  heat  at  higher 
temperatures  where  the  magnitude  of  the  specific  heat  is  large,  no  attempt  has 
been  made  to  fit  the  specific  heat  curve  below  20*K.  The  frequency  distri¬ 
butions  of  the  acoustical  modes  have  been  approximated  by  the  triangular  two- 
dimensional  Debye  distributions,  those  of  the  optical  modes  by  peaked 
Einstein  distributions  or  by  rectangular  distributions  constructed  from  one¬ 
dimensional  Debye  distributions.  In  a  few  cases  the  optical  modes  were 
omitted  or,  more  exactly,  were  assumed  to  have  the  same  frequency  distri¬ 
bution  as  the  acoustical  modes.  It  should  be  noted  that  Tarasov  l12)  has  ap¬ 
proximated  the  specific  heat  in  the  low  temperature  transition  region  by  a 
combination  01  two-  and  three-dimensional  Debye  functions.  Although  it  gives 
a  T2—  T3  transition,  it  is  otherwise  inadequate  in  that  it  fails  to  take  into  ac¬ 
count  the  differences  between  the  various  acoustical  and  optical  modes. 
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2.  LITERATURE  SURVEY 


2.  1.  Method  of  Evaluation 

The  following  procedures  were  used  in  reanalyzing  the  original  data  re¬ 
ported  on  measurements  of  the  specific  heat.  If  the  enthalpy  (total  heat  content) 
was  measured  by,  say,  dropping  a  hot  sample  into  a  cold  calorimeter,  then 
the  experimental  enthalpy  points  were  plotted  against  the  temperature  on 
large  graph  paper.  A  "best  by  eye"  smooth  curve  was  drawn  through  these 
points.  The  specific  heat  at  constant  pressure  was  obtained  from  the  slopes 
of  tangents  to  the  enthalpy  curve  at  a  standard  set  of  temperatures.  Since 
there  is  always  some  error  in  determining  the  slopes  of  tangents  to  a  graphi¬ 
cal  curve,  these  Cp  values  were  plotted  against  the  temperature  and  a  "best 
by  eye"  smooth  curve  drawn  through  them.  This  smooth  curve  was  taken  as 
the  average  specific  heat  of  the  sample,  and  the  reported  tabular  values  are 
the  coordinates  of  this  curve.  In  some  cases  it  was  felt  that  the  scatter  in  the 
enthalpy  points  at  the  higher  temperatures  made  it  impossible  to  adequately 
determine  a  smooth  curve.  Such  regions  were  omitted  and  results  are  given 
here  for  a  smaller  temperature  interval  than  in  the  original  report. 

If  the  specific  heat  was  measured  by,  say,  adding  a  small  amount  of  heat 
to  a  sample  at  temperature,  then  the  experimental  Cp  points  were  plotted 
against  the  temperature.  A  "best  by  eye"  smooth  curve  was  drawn  through 
these  points,  and  the  coordinates  of  this  curve  at  the  standard  set  of  tempera¬ 
tures  are  given  in  the  tables.  Differences  between  the  smooth  curve  reported 
here  and  the  smooth  curve  reported  in  the  original  article  must  be  ascribed 
to  differences  in  taste  as  to  how  to  fit  the  curve.  The  present  author  does  not 
claim  that  his  curve  is  more  accurate  than  that  of  the  original  experimentalist, 
but  the  differences  form  some  indication  of  the  uncertainty  in  the  coordinates 
of  the  smooth  curve. 

At  certain  temperatures,  the  choice  of  which  depended  on  the  nature  of 
the  experimental  curve,  a  rough  estimate  of  the  error  ACp  has  beon  made. 

This  estimate  is  never  less  than  the  originally  reported  error  and  is  greater 
than  the  reported  error  whenever  the  scatter  in  the  experimental  points  was 
so  large  that  the  average  specific  heat  curve  could  not  be  determined  within 
the  original  reported  error.  Although  the  quantity  ACp  was  not  calculated 
from  a  precise  statistical  formula,  it  is  intended  that  it  represents  limits 
such  that  there  is  about  a  50  per  cent  probability  that  the  true  value  lies  within 
the  range  Cp  *  ACp.  The  temperatures  at  which  the  error  was  estimated  were 
chosen  so  that  ACp  should  vary  monotnnically  between  these  temperatures. 

2.2.  Summary  of  Experimental  Measurements  of  the  Specific  Heat  at  Constant 
Pressure  _ 


in  the  following,  individual  samples  are  designated  by  two  sets  of  symbols, 
such  as  CeNG-D55.  The  first  set  (CeNG)  identifies  the  type  of  graphite 
(Ceylon  natural  graphite)  and  the  second  set  (D55)  identifies  the  author  and 
year  of  publication  (DeSorbo,  1955).  The  explanation  of  these  symbols  and 
the  key  to  the  literature  references  are  given  in  Table  1. 
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Table  1 .  Explanation  of  Sample 


Symbol  Type  of  Graphite 

CeNG  •  D55  Ceylon  natural  graphite 
CS  -  D53  National  Carbon  Co.  grade  CS 
Reac  -  B54  A.E.R.E.  reactor  graphite 
AGOT  -  E45  National  Carbon  Co.  grade  AGOT 
(Data  represent  an  average 
of  7  unirradiated  samples) 

Fab  -  J34  National  Carbon  Co.  fabricated 

Fab  -  Nil  Fabricated  graphite 

Fab  -  K1 1  Fabricated  graphite 

CeNG  -  M23  Ceylon  natural  graphite 

Unk  -  W75  Unknown 

CeRt  -  S  24  Ceylon  natural  and  "retort"  graphite 
7087  -  L56  Speer  Carbon  Co.  grade  7087 
GBH  -  L56  National  Carbon  Co.  grade  GBH 

ATJ  -  N60  National  Carbon  Co.  grade  ATJ 

ATJ  -  F60  National  Carbon  Co.  grade  ATJ 

GBE  -  F56  National  Carbon  Co.  grade  GBE 

3474  -  F56  Speer  Carbon  Company  grade  3474 
CaFi  -  W17  Carbon  lamp  filament 
GBH  -  R57  National  Carbon  Co.  grade  GBH 

GBE  -  R57  National  Carbon  Co.  grade  GBE 

3474  -  R57  Speer  Carbon  Co.  grade  3474 

7087  -  R57  Speer  Carbon  Co.  grade  7087 


ignation 


First  Reference 

Author  No. 


De  Sorbo 

13 

II 

14 

Bergenlid 

15 

Estermann 

16 

Jacobs 

17 

Nernst 

18 

Koref 

19 

Magnus 

9 

Weber 

20 

SchlMpfer 

21 

Lucks 

22 

II 

22 

Neel 

23 

Fieldhouse 

24 

II 

25 

II 

25 

Worthing 

26 

Rasor 

27 

II 

27 

II 

27 

1* 

27 

Table  2  gives  the  results  of  the  reanalysis  of  the  experimental  measure¬ 
ments  of  the  specific  heat  at  constant  pressure.  For  brevity  the  estimated 
error  ACp  is  given  as  a  per  cent  of  Cp.  In  two  cases,  only  a  few  measure¬ 
ments  were  made  and  a  smooth  curve  could  not  be  determined.  In  these  cases 
the  measurements  at  temperature  are  reported.  It  is  impossible  to  adequately 
present  all  the  data  of  Table  2  in  a  single  small  graph.  Figures  1  and  2  show 
all  the  measurements  in  the  temperature  ranges  200  to  460°K  and  1000  to 
2000*K. 

The  extent  to  which  the  specific  heat  is  different  for  different  types  of 
graphite  cannot  be  satisfactorily  determined  from  the  present  data.  It  is 
knowrl2*)  that  below  20aK  the  specific  heat  of  natural  single -crystal  graphite 
is  slightly  less  and  that  of  lampblack-based  and  turbostratic  pyrolytic  graphite 
is  slightly  greater  than  that  of  a  good  quality  coke-based  graphite.  De2~'>Kols 
results,'1  shown  in  part  in  Figure  1,  indicate  that  the  specific  heat  of  Ceylon 
natural  graphite  is  as  much  as  9  per  cent  less  than  that  of  artificial  graphite  in 
the  temperature  interval  120  to  300*K.  This  difference  is  over  twice  as  large 
as  the  combined  experimental  errors  of  the  two  curves.  Above  300*K  the 
measurements  on  natural  and  fabricated  graphite  do  not  seem  to  be  signifi¬ 
cantly  different.  Also,  the  difference  between  the  specific  heat  of  different 
grades  of  graphite  seems  to  be  about  the  same  as  the  difference  between  various 
measurements  of  the  same  grade  of  graphite.  This  may  be  seen,  for  example, 
by  comparing  the  results  for  grades  ATJ  and  GBH  (grade  GBH  should  have 
identically  the  same  thermal  properties  as  grade  ATJ)  with  the  results  for  other 
grades.  It  appears  that  within  the  accuracy  of  the  present  measurements  the 
specific  heat  of  all  coke-based  artificial  graphites  is  the  same  ai  all  tempera¬ 
tures  above  a  few  degrees  absolute. 

Table  3  gives  the  average  specific  heat  at  constant  pressure  for  all  coke- 
based  graphites.  Below  300*K  this  curve  is  based  on  the  individual  measure¬ 
ments  of  only  coke-based  graphites.  Above  300*K  it  is  based  on  the  data  for 
all  graphites.  The  average  curve  is  shown  as  a  dotted  line  in  Figures  1  and  2 
and  represents  a  "best  by  eye"  fit  of  the  individual  measurements  on  graphs 
covering  the  entire  temperature  range  with  proper  consideration  given  to  the 
probable  errors  of  the  individual  curves.  The  quantity  AC  given  in  Table  3  is 
an  estimated  error  such  that  most  of  the  reliable  individual  measurements  fall 
within  the  range  Cp  *  AC.  The  subscript  "p"  is  omitted  since  AC  is  also  the 
probable  error  of  the  specific  heat  at  constant  volume  and  of  the  lattice  specific 
heat. 

Table  3  also  gives  the  estimated  values  of  the  lattice  specific  heat  Cj,  of 
the  difference  (Cp-Cv)  between  the  specific  heats  at  constant  pressure  and  at 
constant  volume, of  the  specific  heat  Cc  of  the  electrons,  and  of  the  contribution 
Ct  of  a  thermally  activated  process,  possibly  the  creation  of  vacancies,  which 
occurs  at  very  high  temperatures.  These  quantities  are  important  for  the 
theoretical  interpretation  of  the  specific  heat  of  graphite  and  are  discussed 
further  in  the  following  sections, 

2.  3.  Difference  (Cp-Cv)  of  the  Specific  Heats 

The  difference  between  the  specific  heat  at  constant  pressure  and  the 
specific  heat  at  constant  volume  is  given  by  the  well-known  thermodynamic 
formula 
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Table  2.  Specific  Heat  at  Constant  Pressure 
for  Various  Types  of  Graphite 

(Part  1) 


T 

CeNG  • 

055 

CS  -  053 

Reac  -  U54 

AGOT  -  E45 

Fab  -  J34 

•K 

CP 

ACp 

CP  ACP 

Cp  ACp 

Cp  ACp 

Cp  ACp 

20 

0.0174 

5 

0.0198  3 

0.0178  2 

0.0187  3 

25 

0.0261 

0.0301 

0.0289 

0. 0298 

30 

0.0406 

0.0441  2 

0.0422 

0.045 

35 

0.0575 

4 

0. 0609 

0.0590 

0.063 

40 

0.0765 

0.0795 

0. 0786 

0.082 

45 

0. 0965 

0. 0987 

0.1003 

0.103 

50 

0.119 

0.118 

0.123 

0.127 

55 

0.143 

0.141 

0.147 

0.153 

60 

0.168 

0.165 

0.173 

0.179 

65 

0.195 

2 

0.191 

0.198 

0. 206 

70 

0.223 

0.219 

0. 225 

0.234 

75 

0.252 

0.249 

0.252 

0. 262 

80 

0. 280 

0. 279 

0. 280 

0.291 

85 

0.  309 

0. 309 

0.  30® 

0.321 

90 

0. 339 

0.  341 

0. 338  1 

0. 351 

0.333  1 

95 

0.369 

0.  373 

0.381 

0. 363 

100 

0.400 

1.5 

0.405  1 

0.412 

0.  395 

110 

0.462 

0.472 

0.475  3 

0.463 

120 

0.525 

0.541 

0.535 

130 

0.592 

0.613 

0.609 

140 

0.659 

0.692 

0.686 

ISO 

0.729 

0  ,‘75 

0.766 

160 

0.799 

0.859 

0.849 

170 

0.871 

1 

0.943 

0.934 

180 

0. 945 

1.028 

1.019 

190 

1.021 

1.112 

1.106 

200 

1.097 

1.197  0.7 

1.193 

210 

1.173 

1.282 

1.279 

220 

1.251 

1.367 

1.368 

230 

1.329 

1.453 

1.456 

240 

1.409 

1.539 

1.545 

250 

1.488 

1 

1.625 

1.635 

260 

1.567 

1.712 

1.725 

270 

1.646 

1.800 

1.815 

280 

1.725 

1.887 

1.903 

290 

1.804 

1.973 

1.990 

300 

1.883 

2 

2.060  1 

2.075  1 

Unit*: 

Cp  -  cal/mole  *K; 

ACp  -  par  cent  of  Cp 
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Table  2.  Specific  Heat  at  Constant  Pressure 
for  Various  Types  of  Graphite 

(Part  c) 


T 

•K 


Fab  -  N1 1 
Cp _ AC 


P 


Fab-Kll 

C  AC 
P _ 


P 


82.5 

87.5 
137.9 
231.4 
235.3 


0. 29 
0.32 


Units:  Cp  -  cal/mole  *K; 


0.676 

1.484 

1.5C6 

ACp  -  per  cent  of  Cp 


(Part  3) 


300 

1.83 

350 

2.40 

400 

2.91 

450 

3.  28 

500 

3.58 

550 

3.82 

600 

4.03 

650 

4.23 

700 

4.41 

750 

4.58 

800 

4.74 

850 

4.88 

900 

5.00 

950 

5.09 

1000 

5.17 

1050 

5.  24 

1100 

5.30 

1150 

5.34 

1200 

5.37 

1250 

1300 

1350 

1400 

10  2.11  2 

2.56 

5  2.98 

3.36 

2  3.69 

3.98 

l  4.21  3 

4.58  ? 

4.82  3 

4.89 

4.95 

5.01 

5.06 

1  5.10 
5.14 
5.18 

2  5.22 
5.25 
5.28 
5.31 

5.34  2 


2.09  5 

2.53 

2.91  1 

3.27 

3.58 
3.85 
4.08 

4.28 
4.46 

4.59 
4.71 
4.81 
4.89 
4.97 
5.04 
5.11 
5.17 
5.23 

5.29 
5.35 
5.40 
5.45 

5.50  1 


•  Cp  -cal/molc  *K ; 


ACp  -per  cent  of  Cp 


Units 


g3S8oo8ooo§o§8|8og 
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Table  2.  Specific  Heat  at  Constant  Pressure 
for  Various  Types  of  Graphite 

(Part  5) 


T 

*K 


GBH-R57 
Cn  AC. 


GBE-R57 
C_  AC. 


3474 -R57 
C„  AC_ 


7087 -R57 
Cn  AC. 


1400 

1500 

5.  78 

5 

1600 

5.  89 

1700 

5.98 

1800 

6.06 

1900 

6.13 

2000 

6.19 

2200 

6.28 

2400 

6. 35 

2600 

6.41 

2800 

6.46 

3000 

6.51 

3200 

6.58 

3400 

6.71 

3500 

6.86 

3600 

7.22 

5 

3700 

7.86 

3800 

8.86 

3900 

11.7 

10 

Units: 

Cp  -  cal/molc  * 

5.37 

10 

5.45 

5.53 

5.52 

5.66 

5.58 

5.76 

5.64 

5 

5.83 

5.70 

5.89 

5.75 

5.94 

5.  84 

6.01 

5.92 

6.05 

5.99 

6.10 

6.05 

6.15 

6.12 

6.22 

6.21 

6.32 

6.35 

6.49 

6.50 

6.63 

6. 78 

5 

6.83 

7.31 

7.18 

8.43 

10 

7.80 

9.2 

ACp  -  per  cent  of  Cp 


8.  56  10 

8  5.70 

5.82 

5  5.93 

6.03  5 

6.12 

6.21 

6.34 

6.44 
6.53 
6.60 
6.66 
6.72 
6.81 
6.92 

5  7.10  5 

7.44 

8.06 

9.9  10 
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N-4454 

Figure  1 .  Survey  of  the  Specific  Heat  at  Constant 
Pressure  from  200  to  460*K 
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ie  Specific  Heat  at  Constant  Pressure  from  1000  to  20Q0°K 
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Caiti:  TvmptrtUix*  •  *K.  Sp«cUlc  KtUt  -c*l/mol# 


Cp-Cv=Vp*BT  . 


(li 


where  V  is  the  molar  volume,  p  is  the  volume  coefficient  of  thermal 
expansion,  and  B  is  the  bulk  modulus.  For  a  porous  material  such  as  graphite 
V,  p,  and  B  must  be  evaluated  for  the  crystals  and  not  for  the  bulk,  porous 
material.  All  quantities  must  be  evaluated  at  the  temperature  T. 

2.3.1.  Volume  Coefficient  of  Thermal  Expansion  and  Molar  Volume 

Let  a|,  a2.  a3  be  the  lattice  constants  of  the  unit  cell  of  any  crystal. 
The  linear  coefficient  of  thermal  expansion  is  defined  by 


°i  = 


(2) 


Integration  gives 


T 

ai(T)  =  a1(T0)  exp  [  J  Oj  d  T  ] 
T0 


(3) 


The  volume  of  the  unit  cell  is 


V  =  k  a,  a*  a3 


(4) 


where  k  is  a  numerical  constant.  It  follows  from  (3)  and  (4)  that 

T 


V(T)  =  V(T0)  exp  I  J  (  o,  «•  o2  +  o,)dT]  . 


The  volume  coefficient  of  thermal  expansion  is  defined  by 

M  . 

V  dT  /_ 


On  differentiating  (5)  and  substituting  into  (61  one  obtains 


(5) 


(6) 


P  =ai*u2*aj 


(7) 
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This  result  is  rigorously  correct.  The  often  made  claim  that  it  holds  only  in 
the  approximation  of  neglecting  certain  higher-order  terms  is  an  erroneous 
conclusion  based  on  starting  with  the  formula 


a.  {T)  =  a.(T0)[l+Qi(T-T0)]  , 


(8) 


which  is  itself  only  an  approximation  to  the  correct  formula  (3). 

If  c  denotes  the  interlayer  distance  in  the  graphite  lattice  in  Angstrom 
units  (c  ~  3.  36  A)  and  a  denotes  the  length  of  the  graphite  unit  cell  in  the 
basal  plane  in  Angstrom  units  (  a  ~  2.46  A),  then  the  density  of  the  crystal 
is  given  by 


d  =  {46. 0536/ a*  c )  g/cm3 


(9) 


and  the  crystallite  molar  volume  on  the  physical  scale  of  atomic  weights  is 
given  by 


V  =  (0.260  884  a*c)  cmJ/mole  . 


(10) 


A  brief  literature  survey  has  been  made  of  the  X-ray  values  of  the  lattice 
constants  of  different  types  of  graphite  as  a  function  of  temperature.  The 
linear  coefficients  of  thermal  expansion  were  obtained  from  the  slopes  of  the 
curves  of  lattice  spacing  versus  temperature.  The  c-axis  spacing  and  the 
coefficient  of  thermal  expansion  are  definitely  different  for  significantly 
different  types  of  graphite,  such  as  natural  single-crystal,  coke-based,  and 
lampblack-based  graphites.  The  a-axis  spacing  and  coefficient  of  thermal 
expansion  appear  to  be  the  same  for  different  types  of  graphite.  From  a 
study  of  several  sources average  curves  of  lattice  spacings  and  linear 
coefficients  of  thermal  expansion  have  been  constructed  for  natural  single  - 
crystal  and  annealed  pyrolytic  graphite  and  for  a  good  quality  coke-based 
graphite  with  a  room  temperature  c-spacing  of  3.  360  A.  The  values,  given 
in  Table  4,  have  been  made  self-consistent  in  the  sense  that  a  lattice  spacing 
and  its  coefficient  of  thermal  expansion  satisfy  equation  (3).  The  coefficient 
of  thermal  expansion  curves  are  typical  of  a  class  of  graphites  and  the  values 
for  individual  samples  may  differ  by,  say,  10  per  cent  above  100*K  and  by  a 
factor  of  two  below  100*K.  The  data  above  28Q0*K  were  obtained  by  a  linear 
extrapolation  of  the  coefficient  of  thermal  expansion  curves. 

The  volume  coefficient  of  thermal  expansion  of  graphite  can  be  calculated 
from  the  data  of  Table  4  and  equation  (7)  with  the  condition  that  ai  =  a* ,  and 
the  molar  volume  is  obtained  from  equation  (10). 

2.3.2.  Bulk  Modulus 

The  bulk  modulus  has  been  measured  at  room  temperature  but,  so 
far  as  the  author  knows,  has  not  been  measured  at  higher  temperatures.  An 
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Table  4 


Natural  and  Coke- 
Based  Graphite 
T  a  ai 


0 

2.4619 

0.00 

20 

2.4619 

-  0.01 

40 

2. 4619 

-  0.08 

60 

2.4619 

-  0.15 

80 

2.4619 

-  0.3 

100 

2.4618 

-  0.4 

120 

2.4618 

-  0.6 

150 

2.4617 

-  0.9 

200 

2.4616 

-  1.3 

250 

2.4614 

-  1.5 

300 

2.4613 

-  1.5 

350 

2.4611 

-  1.3 

400 

2.4609 

-  1.1 

450 

2. 4603 

-  0.9 

500 

2.4607 

-  0.6 

600 

2.4606 

-  0.2 

700 

2.  1606 

+  0. 1 

800 

2.4607 

C.4 

900 

2.4608 

0.6 

1000 

2.4609 

0.8 

1200 

2.4614 

1.1 

1400 

2.4620 

1.3 

1600 

2.4627 

1.4 

1800 

2.4634 

1.5 

2000 

2.4641 

1.5 

2200 

2.4649 

1.5 

2400 

2.4656 

>.5 

2600 

2.4663 

1.5 

2800 

2.4671 

1.5 

3000 

2.4678 

1.5 

3200 

2.4686 

1.5 

3400 

2.4693 

1.5 

3600 

2.4700 

1.5 

3800 

2.4708 

1.5 

Units: 

T  -  *K ;  a, 

c  -  A; 

attice  Spacings  and  Crystallite 
aefficients  of  Thermal  Expansion 
Natural  and  Coke-Based  Graphites 


Natural  and 

Coke-Based 

Pyrolytic  Graphite 

Graphite 

c  a  3 

C  Oj 

3.3368 

0.0 

3.3417 

0.0 

3. 3368 

0.3 

3. 3417 

0.6 

3. 3369 

5.3 

3.3419 

6.1 

3. 3375 

12.8 

3.3423 

14.2 

3.3385 

16.7 

3.3437 

17.5 

3.3397 

18.7 

3. 3449 

19.3 

3.3410 

20.1 

3.  3463 

20.7 

3.3431 

21.6 

3.3485 

22.2 

3.3467 

23.4 

3.3523 

24.4 

3. 3507 

24.8 

3.3566 

26.0 

3.3550 

25.9 

3.3610 

26.9 

3.3594 

26.8 

3.  3656 

27.3 

3. 3640 

27.3 

3.3702 

27.6 

3.3686 

27.6 

3. 3749 

27.8 

3. 3732 

27.8 

3.3796 

27.9 

3. 3827 

28.0 

3. 3890 

28.0 

3.3922 

28.2 

3.3989 

28.2 

3.4018 

28.4 

3.4082 

28.4 

3.4115 

?9.6 

3.4179 

28.6 

3.4213 

28.9 

3.4277 

28.9 

3.4413 

29.5 

3.4478 

29.5 

3.4621 

30.5 

3.4684 

30.1 

3.4838 

32.1 

3.4895 

30.6 

3.5070 

34.3 

3.5111 

31.0 

3.5320 

36.7 

3.5331 

31.5 

3.5589 

39.1 

3.5556 

31.9 

3.5877 

41.5 

3.5785 

32.2 

3.6184 

43.9 

3.6017 

32.6 

3.6512 

46.3 

3.6254 

32.9 

3.6861 

48.7 

3.6495 

33.3 

3.7230 

51.1 

3.6740 

33.7 

3.7622 

53.5 

3.6990 

34.0 

3.8036 

55.9 

3,7244 

34.4 

3.8473 

58.3 

3.7502 

34.8 

-  1 0"6  / 

*K 

.  oj 


estimate  of  its  temperature  dependence  can  be  made  in  the  following  way.  The 
bulk  modulus  of  most  materials  varies  approximately  linearly  with  temperature 
except  near  the  absolute  zero,  where  it  is  independent  of  temperature.  Hence, 
it  is  reasonable  to  set 


B(Po,  T)  =  B(Po,  T0 )  + 


(t*t0) 


(id 


where  pQ  and  T0  denote  atmospheric  pressure  and  room  temperature.  The 
temperature  derivative  can  be  estimated  by  a  method  suggested  to  me  by 
Charles  S.  Smith.  If  we  consider  that 


Using  the  definition  of  the  bulk  modulus, 


B=  -V 


(12) 


(13) 


and  the  definition  (6)  of  the  volume  coefficient  of  thermal  expansion,  one  can 
write  (12)  in  the  form 


It  has  been  found  that  in  several  materials  the  first  term  on  the  right-hand  side 
is  at  least  a  factor  of  10  smaller  than  the  second  term;  and  we  will  make  the 
approximation  of  setting  (9B/0T)y  equal  to  zero.  The  pressure  derivative  can 
be  calculated  from  the  compressibility  data. 
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Bridgman  has  measured  the  change  in  volume  of  a  sample  of  Ceylon 
graphite  as  a  function  of  pressure.  However,  after  cycling  the  sample  to 
25,  000  kg/cm2  and  back  to  atmospheric  pressure  the  density  was  only  2. 23 
g/cmJ,  which  corresponds  to  a  porosity  of  1. 64  per  cent.  Presumably  this 
porosity  decreases  as  the  pressure  increases,  thereby  causing  the  measured 
compressibility  to  be  too  large  by, possibly,  10  per  cent.  Kabalkina  and 
Vereshchagin  '*r'  have  measured  the  change  in  c-spacing  of  both  Ceylon  and 
fabricated  graphite  as  a  function  of  pressure.  The  scatter  in  their  data  is 
large  enough  to  mask  any  difference  between  the  two  types. 

The  compressibility  and  its  pressure  derivative  can  be  related  to  the 
change  in  c-spacing  as  follows.  The  linear  and  volume  compressibilities  are 
defined  by 

1  da  \ 

X,  =  '  I  -3p  )  <15> 


1_  dc 
c  dp 


Applying  these  formulas  to  equation  (10),  one  obtains 


X  =  2X|  +  X, 


The  high  pressure  X-ray  data  can  be  represented  by  an  equation  of  the  form 


(c0  -  c)  / Cq  =h,(p-p0)  +  Mp-po)£. 


where  ht  and  h2  are  constants  and  c0  is  the  lattice  spacing  at  atmospheric 
pressure  p0.  Differentiation  with  respect  to  p  yields 

x  s:.*?'  .  „ 

*  1  -h,(p-p0)  -h2  (p-po)1 


At  p  =  p„ 


X5  (P0)  =  h,  , 
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so  hj  is  the  linear  c-axis  compressibility  at  atmospheric  pressure.  Differenti¬ 
ation  of  (20)  with  respect  to  the  pressure  and  evaluation  at  p  =  p0  yields 


(22) 


(23) 


The  X-ray  data  indicate  that  Xj  is  very  small  over  a  large  change  in  pressure, 
so  we  will  make  the  approximation  of  setting  (8Xj  /8 pW  =  0  and  of  using  (22) 
to  give  the  value  of  (8x/8p)y.  Finally,  the  bulk  modulus  by  definition  is  the 
reciprocal  of  the  compressibility, 


B  =  1/x 


(24) 


from  which  it  follows  that 


8B  1  axA 

¥i  :  '?  ~3p~) 

/T  /T 


(25) 


The  volume  compressibility  data  can  be  represented  by  an  equation  of  the 
form 


(VQ-V)/V0  =  g,  (p-p0)  *  g2  (P-Pq)1  ,  (26) 


where  g|  and  g2  are  constants.  An  elementary  but  tedious  analysis  indicates 
that  porosity  causes  a  significant  change  in  g|  and  a  negligible  change  in  g2. 

If  (V0  -  V)/VQ(p  -p0)  and  {  c„  -  c  )  /  c0  { p  -  p0 )  are  plotted  versus  (  p  -p0),  the 
results  should  be  straight  lines  with  intercepts  at  p  -  p0  of  gi  and  ht  and  with 
slopes  of  g2  and  h2.  Bridgman's  smoothed  data  give  an  extremely  linear 
curve  out  to  15,  000  Kg/cm2  and  a  definite  deviation  from  linearity  at  higher 
pressures.  The  scatter  in  the  original  Russian  data  is  too  large  to  confirm 
the  linear  variation.  A  derivation  similar  to  that  leading  to  equation  (22) 
yields 

§£-j  *  gf  *  2g2  at  p  -  p0  .  (27) 

/T 


18 


It  turns  out  that  gf  and  hf  are  about  a  factor  of  10  smaller  than  2g2  and  2h2. 
It  follows  from  this  and  the  approximate  equality  of  the  pressure  derivatives 
of  X  and  Xj  that  g2  and  h2  are  approximately  equal  and,  hence,  the  slopes  of 
the  curves  mentioned  above  should  be  nearly  the  same. 

Since  Bridgman's  data  cover  a  much  larger  pressure  range,  his  results 
were  used  to  determine  g2;  and  we  make  the  approximation  that  h2  =  g2.  A 
"best  by  eye"  line  with  slope  h2  was  drawn  through  the  Russian  data  and  the 
intercept  at  p  =  p0  taken  as  the  value  of  h2.  The  results  are 


X3  =  ht  =  (2.70  ±  0.1)xl0-12  cm2/d 
h,  =  -(32. 71 1  0. 1 )  x  1 0*24  cm4/d2 


(28) 

(29) 


and 


-58. 1  x  10-24  cm4/d2 


(30) 


The  compressibility  Xj  can  be  roughly  estimated  from  the  elastic  constants 
to  be 


X,  =  -0.044  x  10-u  cm2/d  (31) 

with  possibly  a  large  but  unknown  error.  From  (18)  the  volume  compressi¬ 
bility  is  found  to  be 

X  =  (2.61  t  0.  1  )  x  10-‘2  cm2/d  .  (32) 


From  (24)  and  (25)  the  bulk  modulus  and  its  pressure  derivative  are  calculated 
to  be 


and 


B  =  (  3.  83  •  0  1  )  x  10“  d/cm2 


8,5  t  0.6  . 


(33) 


(34) 


From  (11)  and  (14)  one  obtains 

B  =  (  3.  83  -  0.  00077  (T-300)]x  10“  d/cm2  ,  ( !«.) 

in  which  T  is  in  *K.  Within  the  accuracy  of  the  present  measurements  those 
results  apply  to  both  natural  and  fabricated  graphite. 
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2.3.3.  Numerical  Estimate  of  (Cp-Cv) 


The  difference  between  the  specific  heats  at  constant  pressure  and  at 
constant  volume  can  be  calculated  from  equation  (1)  and  the  numerical  data 
given  in  the  previous  sections.  Of  the  factors  in  VpzBT,  only  p  depends 
strongly  on  the  type  of  graphite.  The  values  of  (Cp-Cv)  for  fabricated  graphite 
are  given  in  Table  3.  For  comparison  a  short  set  of  values  for  natural  graph¬ 
ite  is  given  in  Table  5.  Within  a  few  hundred  degrees  of  room  temperature, 
the  values  should  be  correct  to  within,  say,  10  per  cent  but  they  may  be  in 
error  by  a  factor  of  about  2  at  the  higher  temperatures  due  to  a  failure  of 
formula  (11)  to  hold  over  such  a  large  temperature  interval.  Below  1000°K, 
the  difference  between  the  values  for  natural  and  fabricated  graphite  may  not 
be  significant.  At  higher  temperatures,  the  values  for  natural  graphite  are 
definitely  greater  due  to  the  higher  thermal  expansion  of  natural  graphite. 


Table  5.  Average  Values  of  (Cp-Cv) 
for  Natural  Graphite  * 


T 

(Cp-Cv) 

T 

(Cp-Cv) 

0 

0.0 

600 

0.021 

50 

0.00020 

800 

0.030 

100 

0.00161 

1000 

0.040 

150 

0.0029 

1200 

0.050 

200 

0.0043 

1400 

0.060 

250 

0.0058 

1600 

0.073 

300 

0.0076 

2000 

0.  iitij 

350 

0.0098 

2400 

0.147 

400 

0.0120 

3000 

0.205 

450 

0.0141 

3400 

0.237 

500 

0.017 

3800 

0.258 

Units: 

T  -  *K;  (Cp-Cv) 

-  cal/molc  *K 

In  several  studies  of  the  specific  heat  of  graphite,  the  difference  (Cp-Cv) 
has  been  calculated  from  the  formula 

Cp-Cv  =  A  Cp  T,  (36) 

where  A  is  supposedly  a  constant  which  can  be  calculated  from 

A  =  Vpz  B  /  Cpz  (37) 

using  room  temperature  values.  The  constancy  of  A  can  be  checked  by  evalu¬ 
ating  the  right-hand  side  of  (37)  within  a  few  hundred  degrees  of  room  tempera¬ 
ture,  where  equation  (35)  for  B  should  be  reliable.  The  results  are  given  in 
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Table  6  for  fabricated  graphite.  The  complete  lack  of  constancy  of  A  can  be 
traced  to  the  fact  that  Cp  decreases  to  small  values  at  much  higher  tempera¬ 
tures  than  does  p. 


Table  6.  Values  of  Vp2  B  /  Cp2  for 
Fabricated  Graphite 


a 

T 

Vp2  B/Cp2 

50 

411 

300 

6.4 

100 

106 

400 

3.6 

150 

35 

500 

2.6 

200 

16 

600 

2. 1 

250 

9.4 

700 

1.8 

Units: 

T  -  *K;  Vp2  B/Cp2 

-  10“6 

mole/cal 

2.4.  Electronic  Specific  Heat 

The  specific  heat  of  the  conduction  electrons  can  be  calculated  from  the 
electronic  band  structure  of  single -crystal  graphite  under  the  restriction  that 
the  lattice  constants  or,  more  precisely,  the  band  parameters  do  not  change 
with  temperature.  Komatsu  and  Nagamiya*3*'  have  derived  the  basic  formula 
for  the  electronic  specific  heat.  J.  W.  McClure  has  supplied  the  author  with 
a  slightly  more  accurate  version  of  this  formula.  McClure  further  suggests 
that  at  very  high  temperatures  the  electronic  specific  heat  should  approach 
the  values  calculated  for  a  purely  two-dimensional  structure.  Therefore, 
following  McClure  we  take 

Cc  =  3R  (1.  32546)  (kT/yQ) 2  (0.5  +  0.096797  y^kT 

+  0.  501 557  kT/yt )  for  T  <  0.  73394  /  k  (38) 

and 

Cc  =  3R  (1.32546)  (kT/y0) 2  for  T  >  0. 73394  y,  /  k.  (39) 

in  which  R  is  the  gas  constant,  k  is  Boltzmann's  constant,  and  y0  and  y,  arc 
the  band  parameters  in  the  established  notation.  Actually,  formula  (38)  is 
valid  for  T  <  <  y,  /k  and  (39)  is  valid  for  T  >  >  y,  /k;  but  in  the  vicinity  of  T  - 
Y|  /k  the  two  formulas  give  almost  equal  values  and  it  is  sufficiently  accurate 
to  use  each  to  T  =  0.73394  Yl  / k,  at  which  temperature  they  give  equal  values. 
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McClure  suggests  that  the  best  values  for  the  band  parameters  for 
natural  single-crystal  graphite  are 

y0  =  2. 8  eV  and  y,  =  0. 27  eV.  (40) 

With  these  values  (38)  and  (39)  reduce  to 

Ce  =  (2.27  T  +  3.74  x  10*3  T2  +  1.20  x  10"‘  T3) 

x  10~*  cal/mole  *K  for  T  <  2300*K  (41) 

and 

Ce  =  7.48  x  10’’  T 2  cal/mole  *K  for  T  >  2300  #K.  (42) 

in  which  T  is  in  *K.  The  values  of  the  electronic  specific  heat,  computed  from 
these  formulas,  are  given  in  Table  3. 

From  measurements  of  the  specific  heat  below  2*K,  van  der  Hoeven  and 
Keesom'2*'  find  that  the  coefficient  of  the  linear  term  is  3.  30  for  Madagascar 
single-crystal  graphite  and  about  5  for  fabricated  graphites.  It  is  beyond  the 
scope  of  this  work  to  investigate  the  cause  of  the  difference  between  the  values 
2.27  and  3.30,  or  to  investigate  the  electronic  specific  heat  of  the  less -perfectly 
crystalline  fabricated  graphites  at  higher  temperatures.  For  most  temperatures 
the  values  given  in  Table  3  are  sufficiently  accurate.  However,  at  very  high 
temperatures  an  error  of  a  factor  of  two  in  the  electronic  specific  heat  of 
fabricated  graphite  wo>-ld  be  of  interest  in  connection  with  the  possibility  that 
the  lattice  specific  heat  exceeds  the  Dulong  and  Petit  value  in  this  range. 

2.5.  Specific  Heat  Ct  at  Very  High  Temperatures 


On  measuring  the  specific  heat  of  four  grades  of  fabricated  graphite  Rasor 
and  McClelland'2''  found  that  between  3000*K  and  the  sublimation  temperature 
at  about  3920*K  the  specific  heat  doubled  in  value.  Their  data  indicated  a 
thermally  activated  process  with  an  activation  energy  of  7.  7  t  .5  eV/atom. 

By  fitting  their  data  they  arrived  at  the  following  formula  for  the  contribution 
Cg  to  the  specific  heat  of  this  thermal  process: 

Ct  =  5.6  x  10‘7  T'2  exp  (-8. 9  x  104/T)  cal/mole  *K,  (43) 

in  which  T  is  in  *K.  Values  calculated  from  this  equation  arc  given  in  Table  3. 
Although  the  constants  in  (43)  may  not  be  exactly  right  for  the  average  Cp  curve 
used  in  this  report,  the  agreement  is  close  enough  to  confirm  the  general 
correctness  of  the  form  of  equation  (43). 

Rasor  and  McClelland  cite  evidence  both  from  their  specific  heat  and  from 
their  thermal  conductivity  work  that  the  mechanism  of  this  process  is  the 
thermal  creation  of  vacancies.  In  particular  they  note  that  the  activation  energy 
of  7. 7  eV  is  close  to  the  heat  of  sublimation  7. 4  eV,  which  they  take  as  approxi¬ 
mately  the  energy  of  formation  of  a  vacancy.  However,  recent  electron 


microscope  studies  of  vacancy-controlled  processes  in  graphite  by  Baker  and 
Kelly'39'  indicate  that  the  energy  of  formation  of  a  vacancy  is  only  about  3.3  eV; 
so  the  interpretation  of  Rasor  and  McClelland1  s  results  in  terms  of  vacancies 
is  not  completely  certain. 


Lattice  Specific  Heat 


The  specific  heat  at  constant  volume  may  be  considered  to  be  a  function 
of  the  temperature  and  the  volume.  At  atmospheric  pressure  the  volume  of 
the  crystal  is  a  unique  function  V  (T)  of  the  temperature,  so  one  may  write 


Cv  =  Cv  [T.  V  (T)  ]. 


However,  theoretical  calculations  are  usually  based  on  the  concept  of  a  rigid 
lattice  with  lattice  constants  corresponding  to  some  temperature  To,  which 
might  be  room  temperature  or  0°K.  Thus,  Cv  [  T,  V(Tq)]  is  of  greater  theo¬ 
retical  interest,  and  it  is  convenient  to  define  a  quantity  E  (T)  by 


E(T)  =  Cv  [T.  V  (T)  ]  -  Cv  [  T,  V(T0)].  (44) 

The  specific  heat  of  the  lattice  vibrations  of  a  rigid  lattice,  denoted  by 
Cf  [t,  V(T0)  ]  is  even  more  fundamental  and  is  given  by 

C  [T,V(T0)]  =  C  [T,  V(T0)]  +  C  [T,  V(T0)]  +  C  [  T,  V(T0)1 

(45) 

Combining  equations  (1),  (44),  and  (45),  we  obtain  for  the  lattice  specific  heat 
C|  [7,  V(To)]  =Cp[T,  V(T)]  -  VP3BT  -Cc[T,  V(T0)] 

-Ct[T,  V(T0)]  -  E(T).  (46) 

Overton^40)  has  shown  how  E(T)  may  be  calculated  from  experimentally  mea¬ 
sured  quantities.  For  sodium  and  copper,  he  found  that  E(T)  was  smaller 
than  but  not  negligible  compared  to  V|32BT.  Since  the  data  necessary  to  cal¬ 
culate  E(T)  arc  not  available  for  graphite,  this  term  must  be  neglected  here. 


The  values  of  C.  [T,  V(To)]  calculated  by  (46)  arc  given  in  Table  3  for 
fabricated  graphite. 1  The  results  are  not  accurate  enough  for  the  precise 
value  of  Tq  to  matter.  The  lattice  specific  heat  exceeds  the  Dulong  and  Petit 
value  of  5.96  cal/molc  *K  above  about  2300*K.  The  experimental  error  AC 


and  the  approximations  made  in  calculating  C .  from  C  arc  too  large  to  draw 
any  conclusions  concerning  the  existence  of  anharmonic  forces  in  graphite. 


The  values  of  Cf  above  3000*K  arc  not  significant  in  that  Rasor  and  McClelland 
had  to  assume  a  set  of  values  for  Cj  in  this  temperature  range  in  order  to 


derive  their  formula  for  C».  The  computed  values  of  C,  given  in  Table  3  indi¬ 
cate  the  approximate  validity  of  formula  (43)  for  Ct.  Efforts  to  find  analytic 


formulas  to  represent  Cf  [  T,  V(To)j  tire  described  in  the  following  section. 
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3.  ANALYTICAL  REPRESENTATION 
OF  THE  LATTICE  SPECIFIC  HEAT 


3. 1 .  Form  of  the  Frequency  Distributions  in  Graphite 

In  attempting  to  find  an  analytic  function  to  best  fit  an  empirical  curve 
there  are  an  endless  number  of  functions  which  might  be  tried.  To  reduce  the 
work  to  a  reasonable  amount  certain  rather  arbitrary  decisions  have  to  be  made 
at  the  beginning  as  to  the  types  of  functions  to  be  considered.  In  making  these 
decisions  it  is  helpful  to  consider  the  general  nature  of  the  frequency  distribu¬ 
tions  of  the  lattice  vibrations  in  graphite  and  to  review  the  results  of  the  more 
rigorous  Born-von  Kirman  type  calculations. 

There  are  4  atoms  in  the  unit  cell  of  graphite,  so  there  are  12  branches 
of  the  distribution  of  lattice  vibration  frequencies  in  reciprocal  or  wave-vector 
space.  Because  of  lattice  symmetry,  the  branches  are  degenerate  across  the 
top  and  bottom  faces  of  the  Brillouin  zone,  and  it  is  more  convenient  to  con¬ 
sider  that  there  are  6  branches  in  a  Brillouin  zone  which  is  twice  as  high. 

These  6  branches  may  be  classified  into  3  low-frequency  acoustical  modes  and 
3  high-frequency  optical  modes.  Because  of  the  weak  interaction  between 
layers,  the  directions  of  the  atomic  vibrations  are  mostly  cither  parallel  or 
perpendicular  to  the  layer  planes.  This  fact  leads  to  a  further  classification 
of  both  the  acoustical  and  optical  modes  as:  "out-of-plane",  "in-plane  trans¬ 
verse"  in  which  the  vibrations  are  essentially  in  the  layer  plane  and  transverse 
to  the  direction  of  propagation,  and  "in-plane  longitudinal"  in  which  the  vibra¬ 
tions  are  essentially  in  the  layer  planes  and  in  the  plane  defined  by  the  direction 
of  propagation  and  its  projection  on  the  layer  plane  (a  longitudinal  wave  only 
when  the  direction  of  propagation  is  in  the  layer  plane). 

The  restoring  forces  for  the  "out-of-plar.e"  modes  are  weaker  than  for 
the  "in-plane"  modes  and,  correspondingly,  the  frequency  distributions  for  the 
"out-of-plane"  modes  occur  at  lower  frequencies  than  do  those  of  the  "in-plane" 
modes.  Ii  appears  from  the  Born-von  Karman  type  calculations  of  Newell'41' 
Yoshimori  and  Kitano*'42' and  Baldock'41'  that  the  qualitative  shape  of  the 
frequency  distributions  N{v)  is  as  shown  in  Figure  3. 

3.  2.  Representation  in  Terms  of  Debye  and  Einstein  Functions 

From  the  shape  of  the  distribution  functions  as  given  by  Born-von  Karman 
calculations  and  because  of  the  approximate  T2  dependence  of  the  specific  heat 
at  moderately  low  temperatures,  it  appears  to  be  reasonable  to  represent  the 
acoustical  modes  by  two-dimensional  Debye  distributions  for  which  the  specific 
heat  varies  as  T2  at  low  temperatures.  The  optical  modes  should  be  reason¬ 
ably  well  represented  by  rectangular  or,  if  the  distribution  is  narrow,  by 
Einstein  distributions.  The  rectangular  distribution  can  be  constructed  from 
the  difference  of  two  one-dimensional  Debye  distributions.  Although  the  specific 
heat  of  a  one-dimensional  Debye  distribution  varies  as  T  at  low  temperatures, 
the  two  linear  terms  cancel  in  taking  the  difference  of  two  functions,  and  the 
remainder  is  proportional  to  exp  (-0/T)  at  low  temperatures,  as  is  appropri¬ 
ate  for  optical  modes.  It  should  be  noted  that  we  ar-?  using  one-  and  two- 
dimensional  distributions  only  to  approximate  the  true  distributions  and  arc  not 
considering  the  graphite  crystal  as  being  cither  one  or  two  dimensional. 
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Figure  3.  Qualitative  Shape  of  the  Frequency 
Distributions  for  Graphite 


3.2.1.  General  Lattice  Specific  Heat  Function 

Let  Nj  (v)  be  the  frequency  distribution  of  th»  jth  branch  (j  =  1 ,  •  •  • ,  6), 
normalized  such  that 


j*  Nj  (v)  dv  =  l/j  N0, 


{47) 


where  N0  is  Avogadro's  number.  Also,  let  Dj  (x)  and  D£(x)  be  one-  and  two- 
dimensional  Debye  specific  heat  functions  and  E(x)  be  the  Einstein  specific 
heat  function,  where 


x  =  hv/kT,  (48) 

Mathematical  formulas  and  scries  expansions  for  these  functions  arc  given  in 
the  following  section. 

For  the  present  purposes  we  wish  to  consider  three  types  of  distributions 
and  their  specific  heat  functions: 


Linear 


Nj  (v)  =  No  v/  vj2  for  0  <  v  <  vj  (49) 

=  0  otherwise 

(T)=i  RDjIxj)  (50) 

Rectangular 

Nj  (v)  =  No/ 2(1  -dj)  vj  for  djVj  <  v  <  Vj  (51) 

0  <  dj  <  1 

=  0  otherwise 

C|  (T)  =  [R/2(  1  -  dj)]  [D,(Xj)-dj  D,  (djXj)]  (52) 

Einstein 

N;  (v)  =  1  N  6  (v-  vj)  6  =  Dirac  delta  function  (53) 

J  2  J 

Ct  (T)  =  £  RE(xj).  (54) 


A  computer  program  was  prepared  to  compute  the  total  lattice  specific 
heat.  Different  types  of  functions  could  be  chosen  for  the  6  branches  by 
changing  the  input  data.  Fo.'  a  given  choice  of  functions,  the  limiting  frequen¬ 
cies  were  varied  to  obtain  the  best  fit  of  the  experimental  curve.  With  the 
facilities  available,  it  was  impractical  to  vary  more  than  three  frequencies 
during  the  curve-fitting  process.  Therefore,  the  program  was  made  to  depend 
on  only  three  characteristic  frequencies  or,  equivalently,  on  three  character¬ 
istic  or  Debye  temperatures  ,  (Q>z,  and  @ j.  where 

®i  =  h  vj  /  k,  i  =  1,  2,  3.  (55) 

Thus,  the  specific  heat  of  the  jth  branch  depended  on  some  ®  t  (J)»  where 
t( j)  *1,  2,  or  3.  To  allow  flexibility  the  maximum  frequency  oreach  branch 
was  taken  to  be  some  constant  aj  times  one  of  the  three  vj.  The  computer 
program  computed  the  molar  lattice  specific  heat  from  the  formulas 

6 

C,(T)  =  R  £  Fj  (xj),  (56) 

j=l 

where 

xj  =  3^  ®t(jj/T  (57) 
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Fj  (Xj)  =  fj  D2  (Xj)  for  j  =  1,  2,  and  3  (58) 

Fj  (xj}  =  [fj/d-dj)]  [D^x.J-djD.fd.x.)]  +  e.E(x.) 

for  j  =  4,  5,  and  6  (59) 

and  either  f.  =  0  or  e.  s  0  for  j  =  4,  5,  and  6.  The  input  data  to  the  program 
consisted  of  the  f's,  4' s,  d's,  a*s,  t's,  @’s,  and  a  list  of  temperatures  at 
which  the  value  of  Cj  was  desired.  The  input  data  also  contained  the  corres¬ 
ponding  set  of  experimental  values,  which  will  now  be  denoted  by  C|  exP,  from 
Table  3;  and  the  program  computed  the  difference 

Diff  =  Cjj  (T)  -  CfW  (T).  (60) 

The  curve  fitting  was  done  by  means  of  another  program  which  computed 
a  weighted  mean  square  error  MSE  from 

n 

MSE  =  f  I  wj  [ct  (Ti)  -  Or  <Ti)  ]2  .  (61) 

i=  1 

where  n  is  the  number  of  temperatures  used,  wj  is  a  weighting  factor  included 
in  the  input  data,  and  Cj(T)  is  the  value  computed  by  equations  (56)  to  (59). 

The  input  also  included  a  set  of  increments  A  and  the  number  of  times 
A  is  to  be  added  to  ®j.  The  program  computes  MSE  for  the  original 
then  increments  one  of  the  0's  and  recomputes  MSE,  etc.  The  output  gives 
for  each  value  of  @j  a  matrix  array  of  MSE  values  with  the  rows  indexed  by 
the  values  of  and  tne  columns  indexed  by  the  values  of  @2.  The  procedure 
was  to  start  with  a  coaise  net  of  @  values  and  find  the  general  positions  of  the 
absolute  and  relative  minima  in  the  values  of  the  mean  square  error.  By  re¬ 
peating  the  calculations  with  smaller  increments  A  @  the  @  values  giving  the 
minimum  mean  square  error  were  found  to  the  nearest  5  or  10  degrees,  at 
which  point  the  process  was  stopped. 

The  weighting  of  the  mean  square  error  was  done  implicitly  by  the  choice 
of  the  temperatures  used  and  explicitly  by  the  choice  of  the  weighting  factors 
wj.  Ferhaps  the  most  logical  choice  of  the  weighting  factors  is  the  reciprocal 
of  the  square  of  the  experimental  error  A  C, 

w4  =  1  /  [A  C  (Ti)]2.  (62) 

However,  below  300*K  the  experimental  curve  is  known  so  accurately 
that  the  weights  given  by  (62)  are  overwhelmingly  large;  and  in  this  region 
smaller  weights  had  to  be  assigned  on  the  basis  of  how  close  a  fit  is  desired 
in  this  region.  In  order  to  reduce  the  computer  time  to  a  reasonable  amount, 
only  7  temperatures  were  used.  Since  the  primary  objective  is  to  obtain  a 
reasonable  fit  over  the  entire  temperature  range,  no  effort  was 
an  optimum  fit  in  the  T2  region  below  100*K.  The  values  of  T, 
used  in  the  curve  fitting  are  given  in  Table  7, 


made  to  obtain 
-  exp  . 

Cj  ,  and  wj 
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Table  7.  Values  of  T,  Ciexp,  and  Wi 
Used  for  Curve  Fitting 


T 

c/xp 

w« 

1 

100 

0.400 

2500 

200 

1.190 

*>CAA 

300 

2.061 

2500 

500 

3.62 

277 

700 

4.42 

83 

1000 

5.01 

69 

1400 

5.46 

44 

Units: 

T  -  *K;  C/Xp  -  cal /mole 

*K; 

wi  -  mole2  *Kz/cal2 

3.2.2.  Formulas  and  Scries  Expansions  for  Debye  and  Einstein  Functions 

The  Einstein  specific  heat  function  E  (x)  is  defined  by 

E(x)  •=  x2  cX/  (eX  -  1  )2  .  (63) 

For  x  >  0.  3  the  function  E{x)  was  computed  by  (63)  using  the  standard 
computer  subroutine  for  the  exponential  function.  For  x<0.  3  the  function 

X  (x)  =  x/(cx  -  1)  (64) 

was  computed  from  the  expansion 

«0 

x/(ex-l)  =  1  -x/2  ♦  £  (-1)"”*  Bjn„j  x*n/{2n)!  .  (65) 

n  - 1 

where  the  Bln_,  are  the  Bernoulli  numbers 


B,  =  1/6,  B,  =  1/30,  B,  =  1/42,  etc.; 
and  then  E(x)  was  computed  from 

E  (x)  =  X*  ex. 

The  one 'dimensional  Debye  specific  heat  function  is  defined  by 

x 


d,  (xj  =  |  f  i±y.  - 

x  eJ  eY-  1 


cx  -  1 


(66) 


(67) 


(68) 
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and  the  two-dimensional  Debye  specific  heat  function  is  defined  by 


D‘ / 
0 


yzdy 

ey- 1 


2x 

ex-  1 


(69) 


These  functions  were  computed  from  series  expansions  derived  in  the  same 
manner  as  Debye' s(<4)  derivation  of  the  series  expansions  of  the  original 
three-dimensional  specific  heat  function. 


For  x  >  2 


Di 


2 

x 


l 

n  =  ! 


(1  +  nx)  e"nx 


(70) 


x  n=l 


( 2+2n  x+n*x2)e_nx 
T3 


2x 
ex  - 1 


,  (71) 


where 

£{3)  =  1.202056903. 


For  x  <  2 


1}  (x)  «  1 


go 

1 

n=l 


(-l)n_l  {20-5)5^.!  x2n 
(  2  n  +  1  )  .* 


(72) 


Dz  (x)  =  J 


(-l)"'1  (2n-l)B2n.t  xZn 
(n  +  1)  ( 2n )  • 


(73) 


where  the  B2n_j  are  the  Bernoulli  numbers  given  by  (66).  Eight  terms  were 
used  in  each  summation  and  the  break  point  of  x  =  2  was  chosen  so  that  the  two 
series  expansion  for  Dj  and  Dt  would  have  about  the  same  accuracy  at  the  break 
point.  This  procedure  gives  Di  and  Dj  correctly  to  7  or  mors  decimal  places, 
which  is  rather  more  than  adequate  for  the  present  study. 


3. 3.  Numerical  Examples 

From  6  to  9  parameters  aj  and  dj,  depending  on  the  choice  of  E  (x)  or 
D]  (x)  functions  for  the  optical  modes,  had  to  be  assigned  initially  to  determine 
the  general  form  of  a  frequency  distribution.  Then  the  specific  form  of  the 
distribution  was  obtained  by  choosing  the  0  *s  to  minimize  the  mean  square 
error.  The  program  followed  here  was  to  start  with  a  distribution  whose 
general  form  approximated  the  more  exact  frequency  distribution  shown  in 
Figure  3.  Next,  various  simplifications  in  the  general  form  were  made  to  see 
if  acceptable  fits  to  the  experimental  curve  could  be  obtained  with  simpler 
functions.  Only  a  limited  number  of  distributions  could  be  examined  within 
the  contract  period  and  the  distributions  reported  here  are  not  necessarily  the 
best  that  could  be  found.  Nevertheless,  several  distributions  have  been  found 
that  fit  the  experimental  curve  reasonably  well  and  several  points  of  general 
interest  for  the  analytical  representation  of  specific  heat  curves  have  arisen 
from  this  study. 

Distribution  A 


The  distribution  tried  initially  was  similar  to  that  shown  in  Figure  4, 
except  that  the  maximum  frequency  of  the  acoustical  in-plane  transverse  mode 
was  assigned  the  independent  value  vj.  In  order  to  save  space  the  ordinates 
of  the  distributions  shown  in  Figure  4  and  subsequent  figures  have  not  been 
drawn  to  scale,  and  for  convenience  the  abscissa  has  been  labeled  with  the 
characteristic  temperature  instead  of  the  frequency.  It  appeared  from  the 
calculations  that  the  minimum  mean  aquare  error  MSE  for  this  initial  distri¬ 
bution  would  occur  for  vj  or,  equivalently,  greater  than  0.9  (£pi ;  that  is, 
the  acoustical  and  optical  branches  would  overlap.  This  result  was  unexpected 
in  that  Figure  3  indicates  that  <££3  should  be  about  0.  5  .  At  this  point,  it 

was  decided  to  maintain  some  resemblance  to  the  distributions  of  Figure  3  and 
a  gap  was  arbitrarily  set  between  the  in-planc  transverse  distributions.  The 
final  result  is  the  two-parameter  distribution  shown  in  Figure  4  and  denoted 
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Figure  4.  Frequency  Distribution  A. 
MSE  =  0.447  for  the 
Optimum  Temperatures 
Indicated 
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by  Distribution  A.  On  carrying  out  the  curve-fitting  procedure  with  the  tem¬ 
peratures  and  weights  given  in  Table  7,  the  minimum  mean  square  error  was 
found  to  be  MSE  =  0.447  at  ©i  =  2035*K  and  ©2  =  1895°K.  From  the  defini¬ 
tions  (61)  and  (62)  of  MSE  and  w^  it  follows  that  for  a  mean  square  error  of 
unity  the  computed  curve  is  withm  the  limits  AC  of  the  experimental  curve  on 
the  average.  However,  because  the  computed  and  experimental  curves  cross 
each  other  several  times,  the  difference  between  these  curves  may  be  several 
times  AC  in  some  regions,  even  when  MSE  =  1. 

In  order  to  see  how  sensitively  the  optimum  characteristic  temperatures 
depended  on  the  choice  of  temperatures  and  weights  used  to  fit  the  curve. 
Distribution  A  was  refitted  using  14  temperatures  and  larger  weights  which 
decreased  from  15,625  at  100*K  to  100  at  1200eK.  The  minimum  mean  square 
error  with  these  larger  weights  was  MSE  =  0.805  at  @1  =  2045*K  and  ©2  = 
1890*K.  These  values  of  ©)j  and  ©2  arc  in  very  good  agreement  with  those 
found  from  the  smaller  set  of  7  temperatures;  so  the  smaller  set,  given  in 
Table  7,  was  used  for  fitting  all  other  distributions. 

The  lattice  specific  heat  was  computed  from  (56)  using  the  optimum  values 
©1  =  2035*K  and  ©2  =  1895*K  for  a  much  larger  set  of  temperatures.  Table 
8  gives  the  values  of  the  functions  Fj  (x:)  and  Cj(T)  and  of  the  difference  Diff 
as  defined  by  (60).  The  table  was  typea  directly  from  the  output  tape  of  the 
computer;  the  number  following  the  letter  "E"  is  the  power  of  ten  by  which  the 
preceding  decimal  is  to  be  multiplied,  for  example, 

.48378E-05  =  .48378  x  10"*. 

From  the  column  of  Diff  values  it  is  seen  that  the  computed  and  experimental 
curves  cross  each  other  five  times.  Below  200*K  the  magnitude  of  the  differ¬ 
ence  between  the  curves  is  less  than  0.03  cal/mole  *K  but  the  relative  error 
is  up  to  20  per  cent  and  even  larger  below  10*K,  as  expected.  Above  200*K 
the  relative  error  is  about  2  per  cen*  or  less  and  the  computed  curve  is  within 
the  experimental  range  Cjex^  t  AC  except  around  •100*K  where  it  is  within 
the  range  Cj  p  t  2  AC.  Although  the  computed  lattice  specific  heat  for 
Distribution  A  is  not  the  best  that  could  be  found,  it  should  be  a  satisfactory 
approximation  for  most  applications  requiring  an  analytic  specific  heat  function 
for  the  entire  temperature  range. 

Distributions  B,C,  and  D 


The  first  step  in  simplifying  Distribution  A  was  to  assume  that  the  in- 
planc  transverse  modes  have  the  same  frequency  distributions  as  the  in-planc 
longitudinal  modes.  This  approximation  leads  to  Distribution  B  shown  in 
Figure  5.  The  minimum  MSE  is  0.945  for  ©,  =  2175*K  and  ©1  =  2035*K. 
The  fit  to  the  experimental  curve  is  better  at  high  temperatures  but  worse 
near  and  below  room  temperature,  where  the  larger  weighting  factors  cause  a 
greater  mean  square  error  than  for  the  curve  for  Distribution  A. 

Distributions  C  and  D,  shown  in  Figures  6  and  7,  were  tried  to  sec  the 
effect  of  shifting  the  position  of  the  gap  between  the  acoustical  and  optical  out- 
of-plane  modes.  The  minimum  mean  square  errors  are  1.56  and  1.97, 
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Figure  5.  Frequency  Distribution  B.  MSE  =  0.  945 
for  the  Optimum  Temperatures  Indicated 
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Figure  6.  Frequency  Distribution  C.  MSE  =  1.56 
for  the  Optimum  Temperatures  Indicated 
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Figure  7.  Frequency  Distribution  D,  MSE  =1.97 
for  the  Optimum  Temperatures  Indicated 
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respectively,  and  the  fit  to  the  experimental  curve  is  appreciably  worse  below 
500" K  than  for  Distribution  A.  Below  room  temperature,  about  85  per  cent  of 
the  specific  heat  is  due  to  the  acoustical  modes.  Since  the  curve-fitting  pro¬ 
cedure  has  emphasized  this  temperature  range,  it  is  not  surprising  to  find  that 
the  average  characteristic  temperature  of  the  in-plane  acoustical  modes  is 
about  1300  £  100*K  and  that  of  the  out-of-plane  acoustical  mode  is  about  800  ± 
100*K  for  all  four  distributions,  A  through  D.  On  the  other  hand,  the  upper 
characteristic  temperatures  of  the  optical  modes  vary  by  several  hundred 
degrees  and  cannot  be  accurately  determined  from  the  present  experimental 
data. 

Distribution  E 


In  several  theoretical  calculations  the  four  in-plane  modes  have  been  rep¬ 
resented  by  a  single  Debye  function.  Distribution  E,  shown  in  Figure  8,  has 
been  investigated  to  see  the  effect  of  the  approximation  of  neglecting  the  optical 
modes  or,  more  accurately,  of  assuming  that  the  optical  and  acoustical  modes 
have  the  same  frequency  distribution.  Surprisingly,  this  distribution  gave  the 
best  fit  of  all.  The  minimum  mean  square  error  is  only  MSE  =  0.289  for  <3>i  = 
2165*K  and  =  1735*K.  The  good  fit  must  be  considered  as  fortuitous, 
since  there  is  no  theoretical  justification  for  treating  the  optical  modes  in  this 
manner. 

It  has  been  claimed  in  fitting  a  distribution  such  as  Distribution  E  to  the 
experimental  data  that  the  low  temperature  specific  heat  is  due  almost  entirely 
to  the  acoustical  out-of-plane  mode  and  that  its  characteristic  temperature 
(695* K  in  this  case)  can  be  determined  from  the  low  temperature  data  by  neg¬ 
lecting  the  in-planc  modes  entirely.  This  is  not  the  case.  In  spite  of  the  high 
characteristic  temperature  of  2l65,K  of  the  four  in-plane  modes  of  Distribution 
E,  these  modes  contribute  about  one  third  of  the  specific  heat  at  all  tempera¬ 
tures  below  100®K. 

Distribution  F 


The  simplest  frequency  distribution  investigated  here  is  based  on  the 
approximation  of  assuming  that  all  optical  and  acoustical  modes  have  the  same 
form.  This  approximation  results  in  Distribution  F,  shown  in  Figure  9.  The 
minimum  mean  square  error  is  MSE  =  2. 14  at  (0)|  =  2450*K  and  =  1070*K. 
Table  9,  which  is  similar  to  Table  8,  gives  the  computed  lattice  specific  heat, 
Diff,  F|  (x> )  and  F}(xj)  for  several  temperatures  for  Distribution  F.  Tho  fit 
to  the  experimental  curve  at  high  temperatures  is  as  good  as  for  any  distribu¬ 
tion  investigated  but  the  fit  below  150*K  is  worse.  As  is  to  be  expected  the  fit 
below  150*K  using  the  two  two-dimensional  Debye  functions  of  Distribution  F 
is  much  better  than  the  fit  obtained  by  Magnus  using  two  three-dimensional 
Debye  functions.  For  example,  at  40*K  Magnus'  formula  W  gives  C *  =  0.024, 
compared  to  0.055  for  Distribution  F.  The  experimental  value  is  0.079  t  .002 
cal/molc  *K. 

An  upper  characteristic  temperature  for  the  in-planc  modes  of  around 
2500*K  ha t  often  been  quoted  in  the  literature.  Although  a  value  in  this  range 
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Figure  8.  Frequency  Distribution  E.  MSE  =  0.289 
for  the  Optimum  Temperatures  Indicated 


Figure  9.  Frequency  Distribution  F.  MSE  =2.14 

for  the  Optimum  Temperatures  Indicated 


Figure  10.  Frequency  Dist/ibutfon  G.  MSE  =  0.292 
for  the  Optimum  Temperatures  Indicated 
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was  obtained  for  Distribution  F,  all  the  other  distributions  give  values  closer 
to  2000°K.  It  seems  likely  that  very  little  significance  should  be  giver  to  the 
upper  characteristic  temperature  until  the  shape  of  the  optical  moae.  i  as  been 
determined  more  accurately. 

Distribution  G 


The  narrowness  of  the  distributions  of  the  in-plane  optical  modes  of 
Distributions  A  and  B  suggests  that  these  distributions  could  be  approximated 
by  Einstein  distributions  at  a  single  frequency.  This  point  is  illustrated  by 
Distribution  G,  shown  in  Figure  10,  which  has  three  independent  characteristic 
temperatures.  Probably  because  the  minimization  was  carried  out  with  respect 
to  three  @'s,  a  good  fit  was  obtained  for  Distribution  G,  the  minimum  mean 
square  error  being  MSE  =  0.  292  for  the  optimum  temperatures  shown  in 
Figure  10.  The  computed  and  experimental  curves  cross  seven  times  and  the 
fit  is  particularly  good  in  the  room  temperature  region. 


One  point  may  be  noted  which  applies  to  all  distributions.  Since  the  mean 
square  error  is  a  function  of  Uj  and  vg  2,  it  can  be  plotted  as  a  surface  above 
the  (JjCi  ©2-planc.  This  surface  has  a  single  valley  with  rather  steep  sides. 
The  bottom  of  the  valley  has  a  very  gentle  slope  and  usually  has  one  relative 
minimum  and  one  absolute  minimum.  If  and  }J2  are  changed  from  their 
optimum  values  in  such  a  way  as  to  stay  along  the  bottom  of  the  valley,  then 
rather  large  changes  can  be  made  in  £jpi  and  @2  with  only  a  small  increase 
in  the  mean  square  error.  These  points  arc  illustrated  for  Distribution  C  by 
the  values  of  MSE  given  in  Table  10.  For  example,  the  effect  of  a  35*K  change 
in  (2385  to  2350)  can  be  canceled  by  a  change  of  25*K  in  £p2  (1335  to  1360), 
thereby  producJ.ig  a  fit  which  is  essentially  as  good  as  the  optimum.  The  gentle 
slope  of  the  bottom  of  the  valley  seems  to  explain  the  equal  success  of  elaborate 
calculations  based  on  very  different  models  of  atcn.lc  forces  but  which  have  a 
parameter  that  is  arbitrarily  varied  until  the  MSE  is  in  the  bottom  of  the  valley. 
A  good  fit  does  not  necessarily  confirm  the  correctness  of  the  rest  of  the 
calculation. 


Table  10.  Mean  Square  Errors  for 
Frequency  Distribution  C 


©J 

•K 

©2 

*K 

MSE 

Type  of 
Point 

2500 

1250 

1.92 

Botlom 

2350 

1360 

1.58 

Bottom 

2350 

1335 

1.87 

Side 

2385 

1335 

1.56 

Abs.  Min. 

2385 

1360 

1.81 

Side 

-1950 

-1850 

3.3.1 

Saddle 

-1625 

-2900 

-2.0 

Rel.  Min. 
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4.  SUMMARY  AND  CONCLUSIONS 


Several  conclusions  can  be  drawn  from  the  results  of  the  literature  sur¬ 
vey  of  the  experimental  measurements  of  the  specific  heat  of  graphite.  Also, 
the  survey  indicates  the  need  for  additional  experimental  data  in  several  areas. 

For  engineering  purposes,  it  appears  that  the  specific  heat  of  all  coke- 
base  graphites  is  the  same  at  all  temperatures  and  that  the  average  specific 
heat  curve  is  adequately  known  below  about  2000®K.  Above  2QQ0*K  all  the  reli¬ 
able  data  were  taken  on  the  same  experimental  equipment;  and  some  additional 
measurements  are  desirable  in  this  range,  particularly  of  the  rapid  increase 
in  the  specific  heat  above  3000*K. 

For  scientific  purposes,  the  data  are  less  satisfactory.  The  specific  heat 
of  lampblack-based  graphite,  which  has  been  measured  only  below  20*K, 
should  be  measured  at  higher  temperatures.  This  data  would  provide  a  more 
complete  check  on  theories  which  account  for  the  increase  in  the  specific  heat 
of  this  type  of  graphite  in  terms  of  changes  in  the  elastic  stiffness  constants 
c  jj  and  C44.  There  does  not  appear  to  be  any  very  accurate  (t  1  per  cent) 
measurements  of  the  specific  heat  of  any  type  of  graphite  in  the  range  of  a  few 
hundred  degrees  above  room  temperature.  The  apparent  difference  between 
the  specific  heats  of  natural  and  coke -based  graphite  below  room  temperature 
should  be  investigated  above  room  temperature. 

There  is  a  need  for  more  accurate  data  on  the  difference  between  the 
specific  heat  at  constant  pressure  and  at  constant  volume.  At  present  there 
are  no  measurements  of  the  temperature  dependence  of  the  bulk  modulus  of 
graphite  crystals.  Accurate  values  over  a  temperature  interval  of  even  a  few 
hundred  degrees  would  be  valuable  for  checking  the  estimate  of  the  tempera¬ 
ture  dependence  made  in  this  study.  There  does  not  appear  to  be  accurate 
data  on  the  crystallite  coefficient  of  thermal  expansion  along  the  c-axis  of 
coke-based  and  lampblack -based  graphites  in  the  region  of  a  few  hundred 
degrees  below  and  above  room  temperature,  although  some  data  exist  at  higher 
temperatures.  Until  the  Cp-Cv  term  is  accurately  known,  the  theoretically 
interesting  question  of  an  anelastic  component  in  the  lattice  specific  heat  can¬ 
not  be  answered. 

It  was  found  that  combinations  of  Einstein  and  one-  and  two-dimensional 
Debye  functions  could  be  found  which  both  adequately  approximated  the  lattice 
specific  heat  curve  and  roughly  approximated  the  frequency  distributions  for 
the  acoustical  and  optical  modes.  However,  it  was  found  that  an  equally  good 
fit  to  the  experimental  curve  could  be  obtained  by  a  distribution  which  did  nut 
properly  treat  the  optical  in-plane  modes.  This  result  appears  to  be  fortuitous 
and  the  additional  neglect  of  the  optical  out-of-plane  mode  gave  a  poor  fit. 

The  use  of  two-dimensional  Debye  functions  for  the  acoustical  modes  seems 
to  be  definitely  better  than  the  use  of  three-dimensional  functions  at  tempera¬ 
tures  below  150*K,  but  the  large  experimental  uncertainty  at  temperatures 
above  500*K  makes  it  difficult  to  choose  between  two-  and  three-dimensional 
functions  in  this  region.  Rather  different  sets  of  characteristic  temperatures 
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<®1.  @2>  or  (®i  >  ©2 1  ®j)  could  be  found  that  gave  essentially  equal  fits 
to  the  experimental  curve.  Because  of  this  and  because  very  different  models 
of  the  frequency  distributions  gave  nearly  the  same  result,  it  is  not  possible 
to  determine  unique  characteristic  temperatures.  However,  it  appears  that 
the  highest  frequency  for  which  there  is  a  significant  number  of  in-plane 
optical  modes  corresponds  to  a  temperature  closer  to  20Q0*K  than  to  the  often 
used  2500*K.  The  characteristic  temperature  of  the  acoustical  out-of-plane 
mode  represented  by  a  two-dimensional  Debye  function  is  in  the  range  755  ± 
60*K  for  the  models  which  gave  good  approximations  to  the  experimental  curve. 
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